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We study the transverse field Ising model on a kagome and a triangular lattice using high-order 
series expansions about the high- field limit. For the triangular lattice our results confirm a second- 
order quantum phase transition in the 3d XY universality class. Our findings for the kagome lattice 
indicate a notable instance of a disorder by disorder scenario in two dimensions. The latter follows 
from a combined analysis of the elementary gap in the high- and low-field limit which is shown to 
stay finite for all fields h. Furthermore, the lowest one-particle dispersion for the kagome lattice is 
extremely flat acquiring a dispersion only from order eight in the 1/h limit. This behaviour can be 
traced back to the existence of local modes and their breakdown which is understood intuitively via 
the linked cluster expansion. 
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I. INTRODUCTION 

The interplay of geometric frustration and quantum 
fluctuations in magnetic systems provides a fascinating 
field of research. It gives rise to a multitude of quan- 
tum phases manifesting novel phenomena of quantum 
order and disorder, featuring prominent examples such 
as spin-ice, spin-glass and spin-liquid phases. This rich 
behaviour originates from the frustration of microscopi- 
cally competing magnetic interactions. This can lead to 
a huge number of quasi-degenerate ground-state configu- 
rations and a suppression of magnetic order even at zero 
temperature. The macroscopic ground-state manifold of 
simple model systems is extremely sensitive to perturba- 
tions. Arbitrary linear combinations of the degenerate 
states may form novel ground states, yielding the great 
variety of unusual effects some of which are not fully un- 
derstood to date. 

Perhaps the simplest and hence most tractable mod- 
els featuring such behaviour are the fully-frustrated Ising 
systems in a transverse field [ ]. On lattices like the tri- 
angular and the kagome lattice, the frustration results 
from the incompatibility of the antiferromagnetic bonds 
in a single triangular plaquette. These systems exhibit 
a classical ground-state manifold with a macroscopic de- 
generacy which is lifted in the presence of an infinitesimal 
transverse field inducing quantum fluctuations. 

In the limit of high transverse fields all lattices realize a 
non-degenerate polarized phase. The situation for weak 
transverse fields however is much less trivial. A large 
class of these frustrated transverse field Ising systems ex- 
hibits a nontrivial quantum ordered phase arising from 
classical disorder, a scenario which is also known as order 
by disorder [ - ]. In contrast, the kagome lattice plays 
a distinguished role as it seems to be very reluctant to 
order [ - ]. It is believed to be a quantum paramagnet 
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for arbitrary transverse fields providing a novel instance 
of disorder by disorder in which quantum fluctuations 
select a quantum disordered state out of the classically 
disordered ground-state manifold. There is no conclu- 
sive evidence excluding a possible ordering of the kagome 
transverse field Ising model (TFIM) and its global phase 
diagram remains unsettled. So far, the only known Ising 
system featuring such an unusual disorder by disorder 
scenario is the one-dimensional sawtooth chain [o]. As 
a candidate for disorder by disorder in two dimensions 
the kagome TFIM hence is a subject of great theoretical 
interest. 

The central objective of this paper is the analysis of 
the gapped paramagnetic phase for the kagome TFIM to 
check whether it indeed persists for arbitrary field val- 
ues providing an instance of a disorder by disorder sce- 
nario in two dimensions. For this purpose, we investigate 
the energy gap by high-order series expansions about the 
high- field limit. Additionally, we determine the gap of 
the corresponding low-energy mode for an infinitesimal 
transverse field by means of a Lanczos based algorithm 
performed on appropriate finite clusters. Our findings 
suggest that high- and low-field gaps are indeed adiabat- 
ically connected implying that the kagome TFIM realizes 
a disorder by disorder scenario. 

The analoque problem on the triangular lattice is dif- 
ferent. Here, series expansion of the one-particle gap 
about the high-field limit clearly indicates a second-order 
phase transition in the 3d XY universality class between 
the polarized and an ordered phase realized at low fields. 
Our results are in quantitative agreement with quantum 
Monte Carlo simulations [ ]. 

Finally, we show that a major difference between the 
triangular and the kagome TFIM lies in the existence 
of a localized low-energy mode for the kagome lattice 
that is completely flat up to order eight perturbation 
theory. In this order a specific momentum corresponding 
to the \/3 X ^/3 structure is energetically favored. The 
emergence of this local mode and its disappearance at re- 
markably high orders in perturbation theory can be un- 



derstood via the linked cluster expansion. Furthermore, 
we generalize our findings to a larger class of systems. 

This paper is organized as follows. In Sec. II we briefiy 
introduce the TFIM and discuss the physics of the low- 
field and high-field regime. In Sec. Ill we present the 
derivation of the effective one-particle Hamiltonian in the 
high-field limit using the method of perturbative contin- 
uous unitary transformations (pCUTs) which is imple- 
mented via a linked graph expansion. The one-particle 
gap for the triangular and the kagome TFIM is analyzed 
in Sec. IV and in Sec. V. The more technical part dis- 
cussing the existence of flat bands and their breakdown 
is placed in Sec. VI. Finally, Sec. VII provides a conclud- 
ing summary. 

II. MODEL 

We study the two-dimensional Ising model in a trans- 
verse magnetic field on kagome and triangular lattices. 
The corresponding Hamiltonian is given by 

ff= J ^ <Tf<7| + /l^<7f (1) 

<i,j> i 

with Pauli matrices erf acting on site i, the antiferromag- 
netic nearest-neighbour exchange J > 0, and the strength 
of the transverse magnetic field h. 

For the zero-field case = 0, the Hamiltonian de- 
scribes a classical antiferromagnet. The energy of a sin- 
gle bond is minimized by an antiparallel alignment of the 
two neighbouring spins. However, there exists no state 
with all bonds satisfying the antiferromagnetic condition 
due to the frustrated geometry of the lattice. Conse- 
quently, a ground state is given by a configuration where 
every triangle has only two antiferromagnetic bonds of 
which there are infinitely many combinations in the ther- 
modynamic limit. The resulting extensive ground-state 
degeneracy gives rise to a finite entropy density at zero 
temperature [ ] , which is a main characteristic of highly 
frustrated systems. Let us also mention that the classical 
Ising antiferromagnet on the kagome lattice is disordered 
for any temperature. 

Next, consider the introduction of a small transverse 
field h <^ J which induces quantum fiuctuations by fiip- 
ping the z-components of the spins. For any infinitesimal 
field the macroscopic ground-state degeneracy is lifted 
to some degree giving rise to a discontinuity in the en- 
tropy at zero field h = and zero temperature T = 0. 
An interesting question is whether the perturbed system 
continues into a quantum disordered phase (disorder by 
disorder) or whether quantum fiuctuations are selecting 
a quantum ordered state (order by disorder). 

We approach this question for the kagome and the tri- 
angular TFIM by considering the high-field limit {h ^ J) 
which provides an ideal setting for a perturbative ap- 
proach giving also access to the low-field regime. One 
expects a fundamentally different behaviour of the high- 
field gap for an order by disorder or a disorder by dis- 
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FIG. 1: (a) Countour plot of the one-particle dispersion (jj^^{k) 
as a function of the momentum k — (/ci,/c2) on the trian- 
gular lattice for x = 0.1. The red (solid) line corresponds 
to the first Brillouin zone, (b) Illustration of the complex 
phase of the eigen state corresponding to the one-particle 
gap at ^min = (|7r, — Itt). The state exhibits a x y/S- 
structure. Note that the structure for the second gap mo- 
mentum ^min = ~(|^5 ^^^^ ^ V^x v^-structure with 
a different orientation. We stress that the displayed arrows 
do not represent spins. 



order scenario. For the first case the elementary gap 
should close at the quantum critical point separating the 
polarized phase from the ordered phase present at low 
fields. In contrast, disorder by disorder implies the con- 
tinuance of the disordered quantum paramagnet for ar- 
bitrary transverse fields which requires the ground state 
to be adiabatically connected to the low-field limit such 
that the high-field gap never closes. We therefore aim 
at a precise determination of the excitation gap in the 
high-field phase for which we use a high-order series ex- 
pansion. 



III. HIGH-FIELD EXPANSION 

Starting from the high-field limit h J ^ we use a 
pCUT [10, 11] to derive a quasi-particle conserving effec- 
tive Hamiltonian up to high order in perturbation. For 
this purpose, it is expedient to interpret the elementary 
excitations corresponding to single spin flips as quasi- 
particles above the vacuum which is given by the fully 
polarized state. The spin flips are then described in terms 
of hardcore bosons represented by creation and annihi- 
lation operators b- and b]. Besides the usual bosonic 
commutation relations these operators meet the hardcore 
constraint 

bA=blbl=0 (2) 

which means that every site i can be occupied at most 
by one boson as every spin is either 'up' or 'down' (or in 
bosonic language: 'occupied' or 'empty'). 
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The Hamiltonian (1) can then be written as 
§ = Eo + Y.ni + xY,{blbj+ bib] +h.c.) 

= Ho^x{To^ T+2 + T_2) = Ho^xV (3) 

where x = J/2h is the perturbation parameter, 
^0 = —1/2, and fii = bib- is the local density on site i. 
We identify the Ising term as the perturbation V which 
can be decomposed into a sum of operators incre- 
menting (or decrementing) the number of quasi-particles 
by m = {+2,-2,0}. The pCUT transforms the initial 
Hamiltonian (3) into a quasi-particle conserving effective 
Hamiltonian of the form 




m=0 



C(mi, . . .,mn)Trr, 



(4) 



where the second sum runs over all possible 
{mi, m2, . . . , mn} with rrii = {+2,-2,0} and 
fh = ^ • rrii = 0. The latter condition reflects the 
block diagonality of the effective Hamiltonian i^eff in the 
number of elementary excitations allowing to restrict 
the calculations to subspaces with a defined number of 
quasi-particles. The coefficients C(mi,...,mn) have 
been calculated as exact fractions up to high order in 
perturbation [10]. The monomials of operators in 
order n generate virtual fluctuations that involve at 
maximum n different bonds. Due to the cluster addi- 
tivity of the effective Hamiltonian, only linked virtual 
processes (these are processes that comprise bonds on 
simply connected clusters) have a finite contribution. 
Consequently, the matrix elements of the irreducible 
n-particle blocks of the Hamiltonian do possess a linked 
cluster expansion [11-13]. 

We calculated the effective Hamiltonian on the kagome 
(triangular) lattice in the one-quasi-particle (IQP) sector 
up to order n = 13 (n = 11) in the perturbative param- 
eter X. This is done via a full linked-cluster expansion 
taking into accout the 1830 (1784) topologically differ- 
ent graphs. Note that the pCUT provides results for the 
thermodynamic limit that are exact up to the computed 
order. 
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IV. TRIANGULAR LATTICE 

In this section we discuss the one-particle properties 
for the triangular TFIM which is confirmed to realize an 
order by disorder scenario. The triangular lattice can be 
described by a single-site unit cell. As a consequence, 
the IQP block of the effective Hamiltonian is diagonal- 
ized by Fourier transformation giving a single band for 
the one-particle excitations which is fully characterized 
by its dispersion uj^^{k). A contour plot of the resulting 
dispersion for a small coupling x = 0.1 is depicted in 
Fig. 1(a). The minima of the dispersion uj^^ are located 



FIG. 2: (a) Bare series and dlogPade extrapolants of the one- 
particle gap A^^ as a function of x for the triangular TFIM. 
Inset represents a zoom close to the quantum critical point, 
(b) Illustration of the convergence behaviour of the critical 
values Xc resulting from the dlogPade extrapolants [L/M] in 
dependence of the considered order r = L + M + 1. Connected 
points of the same color correspond to extrapolants with con- 
stant L — M. The dashed line represents the quantum Monte 
Carlo result xi^^^^ — 0.303 [ ]. The horizontal dotted lines 
indicate the corresponding margin of error ±0.009. (c) Crit- 
ical exponent u extracted from dlogPade extrapolants shown 
in (b) as a function of r. Dashed horizontal line indicates 
the critical exponent jy^^^^ ^ 0.67 of the 3d XY universality 
class. 
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It is instructive to plot the complex phase of the eigen 
state corresponding to the gap momenta ^min which is 
shown in Fig. 1(b). For the triangular lattice a x 
structure is found and one expects that linear combina- 
tions of such structures are the most preferable ordering 
patterns (see below). 

In order to detect a possible second-order quantum 
phase transition, we use dlogPade techniques to extrapo- 
late the one-particle gap A*^. To this end, various extrap- 
olants [L, M] are constructed where L denotes the order 
of the numerator and M the order of the denominator. 
In Fig. 2 the corresponding results for the one-particle 
gap A*^ are displayed. 

The gap shows a strong tendency to close in a narrow 
window around Xc ^ 0.305. These results are in quan- 
titative agreement with quantum Monte Carlo (QMC) 

simulations giving Xc*^"^^^ = 0.303 ± 0.009 [ ]. The con- 
vergence behaviour of the critical values Xc for different 
extrapolants compared to Xc^^^-' is shown in Fig. 2(b). 
Altogether, our findings suggest a second-order quan- 
tum phase transition in the vicinity of xi^^^"^ = 0.303 and 
Xc 0.305. 

Furthermore, one can extract the critical exponent u 
from the dlogPade extrapolation. The corresponding re- 
sults are shown in Fig. 2(c). Averaging over all dlogPade 
extrapolants [n, m] of order 11 with |n|, |m| > 2 , one ob- 
tains u = 0.708 ± 0.012 which is slightly larger than the 
expected critical exponent zu ~ 0.67 for the 3d XY uni- 
versality class [ ]. Interestingly, if one biases the dlog- 
Pade extrapolants to have a pole at xi^^^"^ = 0.303, the 
corresponding averaged exponent jy^^^^ = 0.658 ± 0.016 is 
much closer to 3d XY. 

The closing of the gap corresponds to a quantum phase 
transition between an disordered paramagnetic phase for 
X < Xc and a quantum ordered phase for x > Xc- The 
symmetry of the ordered phase is reduced according to 
the X ^/3 structure which is consistent with the bond 
ordered phase suggested in Ref. . 

The triangular TFIM therefore realizes an order by dis- 
order scenario. A second-order quantum phase transition 
in the 3d XY universality class is detected by the closing 
of the high- field gap. This is different for the analoque 
problem on the kagome lattice as we will elaborate on in 
the next section. 



V. KAGOME LATTICE 

In this section we aim at determining the one-particle 
gap of the kagome TFIM on the whole parameter axis. 
Since the kagome lattice has a three-site unit cell (see 
Fig. 4(a)), the effective one-particle Hamiltonian H^^^ 
is only reduced to an hermit ian 3x3 matrix 0Jim{ki^ ^2) 
with /,m G {1,2,3} per momentum k using Fourier 
transformation. One therefore has three one-particle 
bands for the kagome TFIM. Due to the lattice sym- 
metries, it is sufficient to calculate just the two matrix 
elements co'ii(/ci, k2) and a;i2(A:i, /C2). The remaining ma- 
trix elements can be obtained by the following relations 









^33(^1,^2) 


= uJii{-k2,ki 


-k2) 


^13(^1,^2) 


= ^12{ki - /C2, 


-k2) 


^23(^1, fe) 


= ^12 (fe - ^1, 


-ki) 









(6) 



There are two dispersive high-energy bands exhibiting a 
finite bandwidth already in first-order perturbation the- 
ory. Most interestingly, the remaining low-energy band 
cjiow is completely flat up to and including order n = 7 



.(7) 
low 



{k) = l-2x^ 2x^ - 3x^ 



243 



2609 



. (7) 



Note that in general the presence of such a fiat one- 
particle band indicates the existence of a local mode. We 
give a detailed discussion on the physics of such local 
modes in Sec. VI. Here we continue with the analysis of 
the one-particle gap. As the low-energy band starts to 
disperse in order n = 8, quantum fluctuations occuring 
in this order lift the degeneracy in momentum-space and 
select a specific gap momentum which is energetically 
favored. The resulting dispersion ^\c>w of order 13 is 
shown for X = 0.1 in Fig. 3. 
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FIG. 3: Lowest one-particle energy band (^io^\k) as a func- 
tion of £ = {ki,k2) for the kagome TFIM at x = 0.1. Red 
solid line indicates the first Brillouin zone. Note the extreme 
smallness of the bandwidth. 
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The locations of the energy gap are ^min = 
±(|7r, — Itt). In order n = 13, one finds the fohowing 
expression for the one-particle gap 
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(8) 



The form of the resulting eigenvector (up to a global 
phase shift) is given by 



1 

71 



(^l,e^i^,e"^i^) 



(9) 



where the three coordinates of v correspond to the three 
sites of the unit cell (see Fig. 4). If one plots the complex 
phase of the full lattice eigen state, one finds a >/3 x >/3- 
structure illustrated in Fig. 4(b). In an analogue fashion 
to the triangular lattice, one would expect that, if the 
kagome TFIM realizes an order by disorder scenario, the 
ordered phase has a >/3 x >/3 structure. Interestingly, 
this is suggested as the ordering pattern in the bond- 
ordered phase of the kagome Ising antiferromagnet in a 
longitudinal field [ ]. 

Furthermore, a high-temperature expansion of the 
classical kagome Heisenberg model displays a striking 
similarity to our findings: Fluctuations (in this case ther- 
mal) up to order n = 7in/3 = l/Tdo not pick any 
specific momentum. It is only in order eight that the 
momentum corresponding to a x >/3-structure is fa- 
vored [15, 16]. Consequently, thermal fiuct nations of the 
classical kagome Ising model behave not unlike quantum 
fluctuations in the kagome TFIM in the high-field limit. 
One might wonder whether the origin of this equivalence 
can be traced back to geometrical relations. This is in- 
deed true as we will discuss in Sec. VI. 

In the following we analyze the one-particle gap on the 
whole field axis in order to clarify whether the kagome 
TFIM realizes a disorder by disorder scenario. At this 
point we assume that binding effects are small for the 
TFIM and, consequently, multi-particle sectors play no 
role for the excitation gap (for any J/h). This assump- 
tion is reasonable since the leading contribution to the 
two-particle interaction is repulsive and of order n = 2 
perturbation theory. 

For the extrapolation it is convenient to perform an 
Euler transformation 



(10) 



1 



which maps the parameter axis x G [0,(X)[ on the finite 
interval u G [0,1]. One expects that the dlogPade ex- 
trapolants in u are more stable within the finite interval 
u G [0, 1] allowing to estimate the gap for x ^ oo which 




FIG. 4: (a) The kagome lattice has a three-site unit cell. We 
define unit vectors ei and 62 as shown on the right-hand side 
of the upper panel. The position of a unit cell is then given 
by (711,722) with Hi G 2. The three sites of the unit cell are 
denoted by {1, 2, 3}. (b) Illustration of the eigen state corre- 
sponding to the one-particle gap. Arrows denote the complex 
phase of this eigenvector. The pattern realizes a x a/3- 
structure. Note that the structure for the second gap mo- 



mentum /cmin = — (|7r, — Itt) is also a X a/3- structure with 



a different orientation. We stress that the displayed arrows 
do not represent spins. 



corresponds to u = 1. The resulting extrapolations of 
the high-field gap are shown in Fig. 6. 

For u < 0.4, all extrapolations as well as the bare se- 
ries are well converged. A closing of the gap in this re- 
gion can therefore be excluded. Note that the bare series 
of A{x) converges very slowly as the absolute value of 
the coefficients grows exponentially with ascending order 
and, additionally, as the sign of the coefficients alter- 
nate. But also far beyond the high-field regime most of 
the extrapolants indicate a clear trend. Though isolated 
extrapolants exhibit instabilities, there is no tendency 
that the gap closes at a certain critical value. More- 
over, most of the extrapolations approach a finite value 
Aoo ~ 0.20 ±0.03 for u = 1. 

From the pCUT result the lowest one-particle mode 
in the high-field limit is characterized by a momen- 
tum ^gap = (Itt,— Itt) whereas the ground state has 

k = (0,0). If the gap stays finite for arbitrary x with 
A{x 00) > 0, corresponding to a disorder by disorder 
scenario, it implies that the ground state and the lowest 
one-particle excitation in the high-field limit can be adi- 
abatically connected to the low-field regime. It would be 
therefore desirable to additionally analyze the low-field 
limit of the kagome TFIM which we focus on next. 
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In the low- field limit, the kagome TFIM is best ex- 
pressed as follows 



J 



where y = h/J. For zero field y = 0^ the Hamiltonian 
describes the classical kagorae Ising antiferromagnet with 
macroscopic ground-state degeneracy. An infinitesimal 
field y <^ 1 lifts this degeneracy leading to an energy 
splitting between states with different momenta k. 

We assume that the two lowest energy levels are char- 
acterized by momenta ^ = (ground state) and ^gap 
(lowest excited state) such that they are adiabatically 
connected to the corresponding states in the high-field 
limit. One can show that the energy gap in the low-field 
limit A{y) (in units of J) and the gap A{x) (in units of 
2h) are related by the following equation 



|A(,)|^^„ = 2JimA(.) 



(12) 



Consequently, the gap A{x oo) can be determined 
by means of first-order degenerate perturbation theory 
about the low-field regime. 

Due to the extensive ground-state degeneracy, one 
has already in order-one degenerate perturbation the- 
ory infinitely many states in the thermodynamic limit. 
We therefore approach this problem by applying first- 
order degenerate perturbation theory on appropriate fi- 
nite clusters of different system sizes in order to obtain 
an estimate of the energy splitting in the thermodynamic 
limit by finite-size scaling. One expects that the energy 
splitting converges to the gap A{x oo) rather rapidly 
with increasing system size because a finite energy gap 
implies a finite correlation length. 

In order to ensure conservation of momentum the clus- 
ters have to exhibit translational invariance which can 
be realized by periodic boundary conditions. Since finite 
systems only allow for discrete values of k depending on 
their sizes and geometrical shapes, the choice of the clus- 
ters plays a crucial role. We are interested in the en- 
ergy having the gap momentum kgap which is realized in 
systems allowing for a x V^-structure. Their corre- 
sponding unit cell contains N = 9 spins. Fig. 5 shows 
the appropriate clusters considered in this paper. 

Technically, we use a Lanczos based algorithm which 
allows to calculate the lowest first-order energy correc- 
tions in subspaces with fixed momentum k. For the case 
of N = 36 spins (ground-state degeneracy g ~ 10^^) we 
need about ul ^ 100 Lanczos steps for an accuracy of 
10-^ 

The obtained energy gaps for the different system sizes 
are displayed in the inset of Fig. 6. Using a 
scaling, one finds for the thermodynamic limit the es- 
timate Ax=oo I'^h ~ 0.18. This value is remarkably con- 
sistent with the value obtained from the high-field expan- 
sion. 
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FIG. 5: The finite-size samples of the kagome lattice with pe- 
riodic boundary conditions which are used for the first-order 
degenerate perturbation theory in the low-field limit. Here 
N denotes the total number of spins. The finite-size samples 
are built by the thick black triangles labeled by an integer 
number. Thin triangles are displayed to visualize periodic 
boundary conditions. 
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FIG. 6: One-particle gap A/2/i as a function of u. Lines cor- 
respond to dlogPade extrapolants in u — while symbols 
represent the first-order degenerate perturbation theory re- 
sults of the low-field limit for the different finite-size samples. 
{Left Inset): Zoom around u = 1. {Right Inset): First-order 
degenerate perturbation theory results of the one-particle gap 
against 1/N . Solid black line represents a linear fit of the data. 
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As a further consistency check, the value of the high- 
field extrapolants for = 1 is biased to the value 
Aoo/2h = 0.18 calculated in the low- field limit. The 
resulting extrapolants are shown in Fig. 7. The en- 
hanced stability of the extrapolation indicates a reason- 
able consistency between the high-field and the low-field 
results. Altogether, the results provide considerable ev- 
idence that the gap stays finite for all values of x. We 
conclude that the kagome TFIM exhibits a disordered 
paramagnetic phase on the whole parameter axis. 
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FIG. 7: Biased dlogPade extrapolations of the one-particle 
gap as a function of u. The extrapolants are biased to 
A(i^=l)/2/i = 0.18 obtained from the 1/A^-scaling in the 
low- field limit. 



VI. FLAT BANDS 

One of the key differences between the kagome and 
the triangular TFIM is the existence of a completely flat 
low-energy band up to order eight perturbation theory 
for the kagome lattice. It is therefore the purpose of this 
section to analyze the physics of such flat bands in more 
detail. This is done first for the kagome TFIM and it is 
later generalized to a larger class of systems. 

In general, the momentum independence of an one- 
particle energy band indicates the existence of a localized 
one-particle mode as any superposition of momentum ba- 
sis states provides an eigen state of H^^^ • A uniform 
superposition in momentum space corresponds to a state 
that is completely localized in the spatial domain. The 
effective Hamiltonian H^^^ of the kagome TFIM there- 
fore gives rise to a local one-particle mode which is an 
exact eigen state of the Hamiltonian up to order 7 in 
X. As the low-energy mode starts to disperse in order 
n = 8, quantum fluctuations occuring in this order lift 
the degeneracy in momentum space and local modes are 



FIG. 8: (a) Localized mode iV^ioc) on a hexagon of the kagome 
lattice. The notion di denotes the sign in front of the one- 
particle states contained in |'0ioc)- Filled dark circles corre- 
spond to the six sites Pij surrounding the central hexagon, (b) 
Example of a destructive interference between the two virtual 
hopping processes (1) (thick dashed bonds) and (2) (thick 
solid bonds) exhibiting the same absolute value t^j'* = t^fj of 
the hopping amplitude. The initial states are given by the 
one-particle states \i) and \i) while the final state of both 
processes is given by Due to the opposite prefactor of \i) 
and \i) in |'0ioc), both processes interfere destructively. 



destroyed in this order for the kagome TFIM. 

The localized mode of the kagome TFIM is given by 
a superposition of one-particle states located on a sin- 
gle hexagon with alternating relative phases as shown in 
Fig. 8(a). According to the displayed numbering a local- 
ized mode on a certain hexagon can be written as 

|^,oc) = ^(|l)-|2) + |3)-|4) + |5)-|6)) (13) 

where \i) denotes a one-particle state with a particle (spin 
flip) located at site i. 

Based on the quasi-particle picture the locality can be 
explained by destructive intereferences of virtual hopping 
processes which prevents the particles from escaping the 
hexagon. 

We start by explicitly discussing the locality of |V^ioc) in 
leading order perturbation theory. The first-order con- 
tribution H^^ of the effective Hamiltonian (4) is given 

by 



(1) 



eff 



xTo 



b 



(14) 



where the sum is taken over all bonds of the lattice 
and the Tq correspond to nearest-neighbour hopping pro- 
cesses on bond b. The state |V^ioc) is then obvisously an 
eigen state of H^^ , since hopping processes escaping the 
hexagon interfere pairwise destructively on each of the 
six sites p*^- surrounding the hexagon (see also Fig. 8a). 
Such destructive interferences are the guiding principle 
to understand the locality of |'^/^ioc) also in higher orders, 
which is best understood via the linked cluster expansion. 

Indeed, due to the linked cluster theorem, the general 
expression of the effective Hamiltonian in order n can be 
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decomposed as a sum of linked virtual hopping processes 
as follows 

^iff^ = ^" E ^('^) E ■ ■ • 
= ^" E E E ^('^) ^^i^-. • • ■ 

Csub c,„b=ur=i{''i}™=o 
= ^" E r£L , (15) 

C'sub 

where the last sum is taken over all connected subclusters 
C'sub built by the collection of linked bonds {bi}. Note 
that the same bond can appear several times in the same 
operator product t^^t^^ . . . where r!^. acts on bond 
bi and rnj G {2,-2,0} refers to a pair creation, pair anni- 
hilation, or single-particle hopping process. The operator 
rl?^ represents all virtual fluctuations of the connected 

<^sub ^ 

subcluster Cgub in order n. 

Here we are interested in the single-particle hopping 
amplitudes Uj describing a hopping from site i to site j 

tij = {j\H,^\i)-S,^j{0\H,f,\0) , (16) 

where |0) denotes the zero-particle state being equivalent 
to the fully polarized state. Consequently, the contribu- 
tion of a connected subclusters Cgub to such a hopping 
element is given by 

= E(OiraLK)-^M(o|raijo)) 

n 
n 

In the following we visualize a contribution tfj'''^ by col- 
oring the bonds bi of the subcluster Csub- Most impor- 
tantly, the tfj""^^ only depend on the topology of the sub- 
cluster Csub- 

The latter property is crucial to pinpoint the properties 
of the local mode |V^ioc)- The dispersion emerging in order 
n = 8 is then related to non-trivial processes that do not 
occur for n < 8. In essence, all processes up to order 
n = 7 cancel out. This can be understood as follows. 

In order to escape the hexagon any particle has to tra- 
verse at least one of the sites p^j enclosing the hexagon. 

Now, envision a process t\^^ on a subcluster (1) as shown 
in Fig. 8(b). The initial state \i) is assumed to have a 
prefactor +1 in |V^ioc) and one considers the particle to 
escape through the site p^-, to the final state A pos- 
sible counter process on a subcluster (2) is then given by 
mirroring only the part of subcluster (1) being "below" 
the site p^-, (one can always perform a global rotation of 
the lattice such that p*j is located on top of the hexagon). 
In this case the topology of the subcluster (2) is equiva- 
lent to the subcluster (1) and, as a consequence, one has 

= t[fj. Since the prefactor in |V^ioc) of the state is 
— 1, the sum of both processes contained in Eq. 15 cancel 
out when acting on the state |V^ioc) (see Fig. 8(b)). A 




FIG. 9: (a) Left figure illustrates a hopping process t-j""^^ 
on a subcluster Csub visualized as thick solid bonds. The 
right part represents the counter process tf,Y^ on Csub- The 
dashed vertical line is the symmetry axis through site p* when 
acting on |'0ioc)- (b) Example of a hopping process that is not 
annihilated by a counter process. Any counter cluster has a 
different topology and therefore a different contribution. This 
leads to a finite hopping amplitude of a particle leaving the 
hexagon. 



second more complex example of such counter clusters is 
shown in Fig. 9(a). 

There are also subclusters that do not possess counter 
clusters of the same topology as exemplified in Fig. 9(b). 
These subclusters contain virtual fluctuations that are 
not completely annihilated by counter processes when 
acting on IV^ioc)- Hence, they correspond to fluctuations 
that lead to a breakdown of the local mode. 

The next crucial step is to identify the lowest pertu- 
bation order for which a hopping process exists which 
has at least one counter cluster with a different topol- 
ogy. In general those graphs are related to closed loops 
comprising two different p*^- that do not lie on the same 
symmetry line as any counter cluster cuts such a loop. 
Clearly, the smallest nontrivial process with a loop is of 
order n = 6 as shown in Fig. 10(a). In this order, the 
hopping process on this subcluster is always local {i = j) 
for the TFIM and hence does not contribute to the dis- 
persion. For n = 7 the only process with a loop, where 
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the final state \j) is not part of |'0ioc) (and tiierefore could 
lead to a dispersion), is displayed in Fig. 10(b). But, this 
process has the special property that the site j is located 
on a symmetry line of the lattice (one finds j = j' for the 
counter cluster). 

In compliance with our results, the first hopping pro- 
cess leading to a breakdown of |V^ioc) is indeed of or- 
der n = 8. Figure 10 depicts one such hopping process 
^_subA ; ^j^^ corresponding Cgub- Besides this subclus- 
ter, there are also other processes leading to a dispersion 
in order n = 8. All of them are related to non- trivial 
processes with loops. 



(a) , (b) 




n = 6 n = 7 



(c) 




n = 8 

FIG. 10: (a) Illustration of the only looped process in or- 
der n = 6 as thick bonds. The particle in the inital and 
final state is located at the same site = Hence, this 
process does not conrtibute to the dispersion, (b) Hopping 
process contributing for orders n > 7. The corresponding 
counter process is represented by thick blue dashed bonds, 
(c) Hopping process contributing for orders n > 8. For this 
case any counter cluster has a different topology leading to 
the breakdown of the local mode |'0ioc)- 

Based on these geometrical observations for the 
kagome TFIM one finds the following criteria for the 
existence and the breakdown of local modes on a given 
lattice. Consider an arrangement of one-particle states 
located inside a restricted area jC. The sites enclosing 
this local area are denoted by p*^, . It is clear that any 
subcluster leaving C has to pass through at least one p*^, 
which one can think of as a gate. In order to allow de- 



structive interference every gate p*^, has to be connected 
to at least two sites belonging to jC. This restricted area 
provides a local mode up to a certain order n, if it allows 
for a superposition of one-particles states with alternat- 
ing sign ±1 inside £, such that all hopping processes 
^Csub,(m) subclusters Cgub for m < n have a topo- 

logically equivalent counter cluster. For the triangular 
lattice, it turns out that there is no local arrangment 
that complies to the stated criteria. 

General construction rules for localized magnon states 
are already well known from the intensive study in the 
context of macroscopic magnetization jumps [ ]. Our 
geometrical observations are fully consistent with the 
general rules stated in these works. Beyond these con- 
struction rules, we are able to identify single processes 
that lead to a dispersion which also allows to determine 
the exact order of locality on the basis of a quasi-particle 
picture. The dispersion of a localized mode is affected by 
virtual processes that close a loop between different p*^ . 
The p*j are sites enclosing the local mode located in jC. 
Two sites p*^- are considered as different if they do not 
share the same sites in C. 

Let us illustrate these principles by means of another 
TFIM on a frustrated lattice. For instance, consider the 
checkerboard lattice shown in Fig. 11. Here the local- 
ized state IV^foc) is given by an arrangement of particles 
with alternating phase located at the sites of a square 
plaquette. The displayed process has a non- trivial loop 
containing two different sites p^j leading to a breakdown 

of IV^foc)- ^^^^ TFIM the corresponding 

hopping process first appears in order n = 5. Calcula- 
tions for the checkerboard TFIM indeed show that the 
dispersion of the lowest one-particle mode is completely 
fiat up to and including order n = 4 in accordance with 
our predictions. 




FIG. 11: Localized mode |V^f^^J = | (|1) - |2) + |3) - |4)) on 
the square plaquette of the checkerboard lattice. The dis- 
played hopping process from site 1 to site j involving 5 thick 
solid bonds leads to a dispersion as it closes a loop between 
the two different sites pli and pl2- 
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VII. SUMMARY 

We have studied the triangular and the kagome TFIM 
at zero temperature. Such highly- frustrated Ising models 
represent the simplest examples displaying an intriguing 
interplay of classical frustration and quantum fluctua- 
tions. Indeed, the tranverse magnetic field represents a 
singular perturbation on the extensive number of ground 
states present for the classical Ising problem. 

Generically, one expects an order by disorder sce- 
nario for such problems [ , ]. The corresponding zero- 
temperature phase diagram then consists of two phases, 
an ordered one at low fields and a polarized phase at 
high fields separated by a quantum phase transition. Al- 
ternatively, a fascinating disorder by disorder scenario is 
possible. Here no symmetry-broken phase is realized and 
a polarized phase is present on the whole parameter axis. 

In this work we give strong evidence that the kagome 
TFIM displays a notable first instance of an disorder by 
disorder scenario in two dimensions. In contrast, the tri- 
angular TFIM falls in the more conventional class of sys- 
tems showing an order by disorder. 

Our findings are based essentially on high-order se- 
ries expansions of the one-particle gap about the high- 
field limit. To this end a full linked graph expansion 
is performed using the pCUT method. For the trian- 
gular TFIM extrapolations of the one-particle gap give 
convincing evidence for a gap closing and therefore for 
a second-order quantum phase transition. The location 
of the critical point is in quantitative agreement with 
quantum Monte Carlo simulations [ ]. Furthermore, the 
obtained critical exponent is fully compatible with the 
suggested 3d XY universality class. 

In contrast, the extrapolation of the one-particle gap 
for the kagome TFIM displays no tendency to close. Most 
of the extrapolations suggest a finite gap of the order 
« 0.2 in units of 2h even in the low-field limit, which is 
consistent with first-order degenerate perturbation the- 
ory about the classical Ising antiferromagnet. The latter 
is done for finite systems yielding a gap ~ 0.18 in units 
of 2h in the thermodynamic limit. Put together, we have 
determined the one-particle gap on the whole parameter 



axis. Strong evidence is found that the low-field and high- 
field limit are adiabatically connected. Consequently, the 
kagome TFIM realizes a disorder by disorder scenario. 

Let us also address similarities and differences between 
the impact of quantum fluctuations (in this work induced 
by the transverse field) and of thermal fiuctuations on 
classical magnets on the kagome lattice. Whereas ther- 
mal fluctuations also do not induce order for Ising spins, 
the situation is different for the Heisenberg case. Here 
the high-temperature expansion has an analogous struc- 
ture, in particular with regard to the momentum space 
degeneracy up to order eight in J/T [ ]. There, up to 
Mermin- Wagner physics, dipolar, quadrupolar, and oc- 
tupolar orders appear to be induced by thermal fluctua- 
tions [15, 16, 18-22]. 

The momentum independence of the low-energy band 
in the flrst seven orders of perturbation theory for the 
kagome TFIM can be traced back to the existence of local 
modes in the spatial domain. It is only in order eight 
that quantum fluctuations occur which do not interfere 
desctructively and therefore lead to the breakdown of 
the local modes. We have identifled a mechanism for the 
(dis) appearance of this degeneracy which readily applies 
in a generalised setting. 

Our flndings might inspire future research about sim- 
ilar questions in other classes of problems beyond trans- 
verse fleld Ising models. One fascinating idea would be to 
establish a link between the importance of flat bands in 
magnetic systems discussed in this work and the physics 
of fractional Chern insulators also based on the existence 
of flat bands on the one-particle level. 
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